Search for: All records

Many astrophysical applications require efficient yet reliable forecasts of stellar evolution tracks. One example is population synthesis, which generates forward predictions of models for comparison with observations. The majority of state-of-the-art rapid population synthesis methods are based on analytic fitting formulae to stellar evolution tracks that are computationally cheap to sample statistically over a continuous parameter range. The computational costs of running detailed stellar evolution codes, such as MESA, over wide and densely sampled parameter grids are prohibitive, while stellar-age based interpolation in-between sparsely sampled grid points leads to intolerably large systematic prediction errors. In this work, we provide two solutions for automated interpolation methods that offer satisfactory trade-off points between cost-efficiency and accuracy. We construct a timescale-adapted evolutionary coordinate and use it in a two-step interpolation scheme that traces the evolution of stars from zero age main sequence all the way to the end of core helium burning while covering a mass range from 0.65 to 300M_⊙. The feedforward neural network regression model (first solution) that we train to predict stellar surface variables can make millions of predictions, sufficiently accurate over the entire parameter space, within tens of seconds on a 4-core CPU. The hierarchical nearest-neighbor interpolation algorithm (second solution) that we hard-code to the same end achieves even higher predictive accuracy, the same algorithm remains applicable to all stellar variables evolved over time, but it is two orders of magnitude slower. Our methodological framework is demonstrated to work on the MESA ISOCHRONES ANDSTELLARTRACKS(Choi et al. 2016) data set, but is independent of the input stellar catalog. Finally, we discuss the prospective applications of these methods and provide guidelines for generalizing them to higher dimensional parameter spaces. 

ABSTRACT Near-infrared (NIR) observations of normal Type Ia supernovae (SNe Ia) obtained between 150 and 500 d past maximum light reveal the existence of an extended plateau. Here, we present observations of the underluminous, 1991bg-like SN 2021qvv. Early, ground-based optical and NIR observations show that SN 2021qvv is similar to SN 2006mr, making it one of the dimmest, fastest evolving 1991bg-like SNe to date. Late-time (170–250 d) Hubble Space Telescope observations of SN 2021qvv reveal no sign of a plateau. An extrapolation of these observations backwards to earlier-phase NIR observations of SN 2006mr suggests the complete absence of an NIR plateau, at least out to 250 d. This absence may be due to a higher ionization state of the ejecta, as predicted by certain sub-Chandrasekhar-mass detonation models, or to the lower temperatures of the ejecta of 1991bg-like SNe, relative to normal SNe Ia, which might preclude their becoming fluorescent and shifting ultraviolet light into the NIR. This suggestion can be tested by acquiring NIR imaging of a sample of 1991bg-like SNe that covers the entire range from slowly evolving to fast-evolving events (0.2 ≲ sBV ≲ 0.6). A detection of the NIR plateau in slower evolving, hotter 1991bg-like SNe would provide further evidence that these SNe exist along a continuum with normal SNe Ia. Theoretical progenitor and explosion scenarios would then have to match the observed properties of both SN Ia subtypes. 

ABSTRACT The progenitor systems and explosion mechanism of Type Ia supernovae are still unknown. Currently favoured progenitors include double-degenerate systems consisting of two carbon-oxygen white dwarfs with thin helium shells. In the double-detonation scenario, violent accretion leads to a helium detonation on the more massive primary white dwarf that turns into a carbon detonation in its core and explodes it. We investigate the fate of the secondary white dwarf, focusing on changes of the ejecta and observables of the explosion if the secondary explodes as well rather than survives. We simulate a binary system of a $$1.05\, \mathrm{M_\odot }$$ and a $$0.7\, \mathrm{M_\odot }$$ carbon-oxygen white dwarf with $$0.03\, \mathrm{M_\odot }$$ helium shells each. We follow the system self-consistently from inspiral to ignition, through the explosion, to synthetic observables. We confirm that the primary white dwarf explodes self-consistently. The helium detonation around the secondary white dwarf, however, fails to ignite a carbon detonation. We restart the simulation igniting the carbon detonation in the secondary white dwarf by hand and compare the ejecta and observables of both explosions. We find that the outer ejecta at $$v~\gt ~15\, 000$$ km s−1 are indistinguishable. Light curves and spectra are very similar until $$\sim ~40 \ \mathrm{d}$$ after explosion and the ejecta are much more spherical than violent merger models. The inner ejecta differ significantly slowing down the decline rate of the bolometric light curve after maximum of the model with a secondary explosion by ∼20 per cent. We expect future synthetic 3D nebular spectra to confirm or rule out either model.